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ABSTBACT 


Constructive existence and uniqueness results for boundary value problems 
associated uith some simple special cases of the second order equation 
'i** - f(K.y.y'')» 0 < X < are sought. The approach ue consider is to 
convert the differential equation and boundary conditions to an integral 
equation via Green's functions, and then to apply fixed point and contraction 
map principles to a sequence of successive approximations. The approach is 
tested on several applied problems. Difficulties in trying to prove general 
theorems are discussed. 

1. INTRODUCTION. (General Problem of Interest.) 

The topic of interest is the establishment of constructive existence and 
uniqueness for nonlinear, tuo-point boundary value problems 

(1.1a) y’' - f(x,y,y'), 0 < x < 1, 

(1.1b) aiy(0) ♦ a 2 y'( 0 ) = Ci, biy(l) + b 2 y'(l) * C 2 . 

By constructive existence and uniqueness, ue mean results whose proofs suggest 
a method for computing the solution numerically. The type of approach considered 
involves one of many already known methods for obtaining solutions numerically. 
The equation and boundary conditions are converted to an integral equation via 
Green’s functions. Then the solution of the integral equation is sought using 
successive approximations. However, the equation is first parametrized in a 
special way which depends on conditions assumed to be satisfied by f. Some new 
existence and uniqueness results can then be obtained. 

For the initial value problem associated uith (la), there exists theory that 
assures a unique local solution for a large class of equations. However, the 
question of existence and uniqueness or just existence for boundary value 


problras assooiatsd uith (!•) is hard to antuor uniats vary strong assuaptions 
are made on f. Even in aany raoant papers on existenoe and uniqueness t 
assumptions made on f usuall include a oonditi<m (tuoh as a Lipsehita condition* 
a boundedness eonditxon* or nondeoreasing behavior) in the argissant y for all y. 
Thus* although they may contain more general results with ueaker conditions on f 
than in earlier comprehensive works (such at Keller tll> and Bailey* Shampine 
and Ualtman t2l) they can still be too strict for sNuiy applications. Also* 
proofs are not always constructive in nature* and in many practioal applications 
solutions are computed without establishing axistence. 

Our basic idea is to seek constructive existence and uniqueness for various 
eases of (1) bated on assumptions on f that hold only for y(x) that satisfy a 
certain constraint. Although our primary objective is to have results which 
enable t!ie solution to be computed numerically based on the constructive nature 
of the proofs* we have a secondary goal. Ue would like numerical methods to 
lend themselves to the application of results from research in vector acceleration 
of sequences and series. 

2. MOTIVATIOH. 

In the interest of applied mathematics* ue can show problems which arise 
naturally in the applied sciences and which exemplify the general problem (1.1). 
Examples of problems (1.1) can be shown to arise for instance from problems in 
heat transfer* problems in the analysis of chemical reactions* and problems in 
flttid mechanics. In particular* ue mention the following. 

In tlie analysis of the stagnation point shock layer (31* it is shown that 
the total cnthalt.’ is governed by 

(2.1a) y''‘'(x) * KR(y(x))*' - xRy^(x)* nil* 0 S x S 1* 

(2.1b) y(0) = 0* y(l) » 1. 

The quantities K and R are positive constants characterized as a radiation loss 
parameter and a Reynolds number respectively. 

In a problem concerning the analysis of heat and mass transfer in a porous 
catalyst (4)> the following boundary value problem is obtained. 

/ Yfl(l-y) \ 

(2.2a) y-'' = Of y exp j ' 

I 1 + p(l-y) 
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(2.2b) y'(P) • y(l) ■ 1. 

The quantities y, 0 and « are positive constants representative of 
diaiens ionless energy of activation* ,heat evolution and Thiele's Modulus 
respectively. 

Ue also consider vector oases of problee (1.1). In particular* for the 2 
by 2 case* y(x) > (yi(x)* ys(x))^ * and f * (fi* fx>^ • 2 general fere 

of the linenr separated boundary conditions uould be 

Aiy(0) * kzy'tO) - Cl* Bty(l) * Biy'^(l) > Ct 

T 

uhere At* Bi and Bj are two by two matrices and Ci = (cn* cti) and Ct * 
(cit> ott)^ • Examples can be shoun to arise in fluid mechanics. For instance* 
consider the following. 

The unsteady squeezing of a viscous fluid between two parallel plates is 
discussed in (51. With the normal velocity prescribed* the unsteady Navier- 
Stokes equations admit a similarity solution. The similarity equation for the 
axisymmetric case is 

(2.Sa) S(xf'^' ♦ 3f'^ - ff''^) * i"", 

(2.3b) f(0) * f'^(O) » 0* f(l) * I* f'(l) » 0. 

which can be written as a vector case of ( 1 . 1 ) > 

\ _ / V2 

Vz" I ^ S(xyt' ♦ 3yi - yivt') 

(2.^b) yi(o) = 0* ytO) = 0* yi(i) ■ i.-yi^’d) - o. 

where yi = f and yt = f". 

For more examples* see Appendix A. 

3. A BASIC RESULT. 


y 


(2,<ia) 


Investigation into the construction of theorems for various subclasses of . 
problem (1.1) is based on a result in 1 6 ) established for the problem 
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(3.1a) 


y" ■ {(K>y) 


0 S X S 1» 


(3.1b) y(0) » 0, y(i) ■ 0. ORfGffML PAGE 13 
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uhara the gradiant y^ doas not appaar axplioitly. This rasult is an axtansion 
of a siailar rasult dua to Kallar III. First, subtract k>y fron both sidas of 
(3.1a) and eonsidar tha aquivalant problaa 

(3.2) y" - k*y « f(x,y) - k*y, y(0) ■ y(l) ■ •. 

Then, for k> # 0, (3.2) can ba eonvartad by tha Groan's function procadura 

into tha equivalent integral aquation 

(3.3a) y(K) « r gk(K,5) (k* y(?) - f(|,y(?))) d| 

uhara 

1 sinh kx sinh k(l-(), 0 i x < If, 

(3.3b) gk(x,|) ® — 

k sinh k sinh k(l'-x) sinh k|, I < x i 1. 

The theorem in (6) that ue are refering to is the follouing. 

Theorem 0. In the boundary value problem (3.1), let df/dy ba continuous 
for all x € (0,11 and all y . Suppose that there exists N > 0 and 820 

such that 0 < 8 < df/dy < N for all x € 10,1] and all y. Then a unique 

solution of (3.1) exists. For k* = (1/2) (8 + N), it is given by the limit 
of the convergent sequence of functions 

(3.4a) y®(x) = 0, 

(3.<ib) y»*’(x) = \ gk(K,?) lk*y"(5) - f(?,y"(5))l d?, m = o,i,***. 

Proof . Let 


(3.5) e"**’(x) * y«**nx) - y»(x) 
and 

(3.6) lie""' II * max |e""'(x)|, m ■ 0, 

0<x<l 

Then for m « 1.2,***, ue can apply the mean-value thaoram to f(x,« :))- 

f (x.y'"* ' (x) ) to obtain 
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t* af 

(S.7) > \ 9m(n.|) (k* (|*y*C|)-«(|)«*(|))l •»(!) 

ay 
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(3.8) k> - (1/2) (S * N) 


and the bounds on af/ay» the bracketed tern in the integrand of (3.7) satisfias 
0 1 |k> - af/dy| 3 (1/2) (N - 8). Therefore, from (3.7). 


1 

|e«*i(K)| i - (N - 8) \ gH(K.f) d{* ||e»|| , 

2 *9 


(3.9) 


1 

» (l/2)(n - 8) - (1 
k* 


cosh k((l/2)-x) 
cosh(k/2) 


) lle»!| 


N - 8 

< (1 

N ^ 8 


1 

) 

cosh(k/2) 


lle-ll , 


1 . 2 . 


Since this relation holds for all k e 10.11. 
(3.10) ^ Sk l|e"*|| 


(3.11) Hu * (1 ). 

N + 8 cosh(k/2) 

Observe that hh < I, and ||e'"^^|| ^ iik"*l|ei|| . Thus {y"*} is a Cauchy sequanea 
in the space of continuous functions on 1 0.1) uith the norm defined by (3.0). 
Therefore, a continuous limit y(K) exists, to which <y"*(K)} converges uniformly. 
Since the order of the limit operation and the integration can be interchanged 
the limit function satisfies the integral equation (3.3). To establish 
uniqueness, let yi(x) and y 2 (K) be tuo solutions to (3.1). Then they both satisfy 
(3.3) for k* = (1/2) (8 * N). By the same analysis that leads to (3.10). 

Ilyi " y2ll ^ kk llyi - yill • since hu < l. ||yi - yill « 0. or yi « y*. 


0. PROGRESS. (Some Neu Results) 
In 171, tuo theorems uere constructed for 
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(4.1a) 

(4.1b) 


■ f(N>y)t 


e $ M s 1, 


y(0) « yot 


y(l) ■ yv 
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Ona thaorea aatablithaa aNiatanea and unlquanaas aaong all y for uhioh |y(K)| i 
aaK (|y«|> tyi|)» x c (0»11. Tha othar astabliahaa axiatanca and unlquanaaa 
among all y for uhieh 0 i y(x) S H uhara yo and yi ara aaauaad nonnagatlva and 
aax (ye* yi> ^ M. Tha eonditiona aaauaad on f diffar alightly. Baoantly (8lf 
ua have obtained raaulta for a more general aubolaaa of (1.1) » namely 


(4.2a) y^' » fi(x,y) ♦ p(x)y^ 

Theorems have been establiahed for three different seta of boundary conditions. 


(4.2b-l) 

y(0) = yo. 

yd) « yi» 

(4.2b-2) 

y'(0) = 0, 

yd) = yi. 

(4.2b-3) 

y(0) = yo» 

y'd) » 0. 


Since the results for (4.2a) together uith (4.2b-l) have as special cases the 
results reported in (7) for (4.1), ue shall shou a theorem and proof for (4.2a) 
together uith (4.2b-l). 

The approach taken on (4.2a) is as follows. Assuming p(x) has a continuous 
derivative, multiply both sides of (4.2a) by the integrating factor where 


1 

(4.3) v(x) = - \ p(?) d5 . 

2 ->0 


Then write (4.2a) as 

(4,4a) u" = F(k,u), 

where 

(4.5a) u(k) = y(x), 

(4.5b) F(x,u(x)) = e*''‘>‘*fi(x,e''*>‘*u(x)) ♦ q(x)u(x), 

(4.5c) q(x) = (p(x)/2)^ - p'(x)/2 . 
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(4.4b) u(0) • u(l) • 

Nou u* can follou an approach siadlar to tha ona for (5.1). First roplaoa 
(4.4a) ulth tha aquivalant aquation 

(4.4) u" - k*u ■ F(x«u) - k*u. 

Than for k* # 0» (4.4) togathar uith (4.4b) can ba convartad into an Intagral 
aquation by tha Graan’s function procadura for tha oparator (d^/dx* - k*)* i.a.y 


(4.7a) 


5 ! 


u(x) « h(x) \ gK(X 2 |) (k2u(|) - r(|yu(())) d|» 


uhere 


(4.7b-l) 


Gktx.O 


1 f sinh k(l-() sinh kx> 

k sinh k ^ sinh k(l-x) sinh k(» 


0 < X < f, 

1 < X i 1> 


and 


(4.7C-1) 


h(x) 


e*vio» sinh k(l-x) ♦ yi sinh kx 

sinh k 


In addition to assuming that p(x) has a continuous derivativa on (0»11> ua 
shall also assume q(x) of (4.5c) is nonnegativa on (0>11. Then* if ua define 

(4.8) 6t “ min q(x)> N 2 » max q(x)* 

0<x<l OixSl 

ue Mill have Sj > 0* and Nj > 0. (Note. Tha assumption q(x) i 0 can ba uaakanad.) 

Theorem lA. In tha boundary value problem uhich consists of aquation (4.2a) 
together uith (4.2b-D* let max (|yel*|ytl> ^ H. Suppose there exists an Ni > 4 
and a St ^ 4 such that 0 < St < dft/dy < Ht for all x c (0.1) and all y such that 
|y(H)| < e''* '*>max{e'''* ®‘,e’''‘ ”}M. x < t0,l|. Suppose further that 0 S ft(x*y) 

< (Ni^St)y. y > 0* and (Ni«Si)y < ft(x.y) < 0* y i 0* for all x c (0*1) and 
all y such that |y(x)| i e'^' ’•*maxf »e*''‘ ' *)M. x < 10*1). Than there exists 
a unique solution of the problem satisfying |yCx)| < a''< *>max(a'''‘ •* ,a*''' ’ •>H* 

X c 10*1 i. Let S = (SftSt} and N » (Nf»Nz) uhere N 2 and Sj are given by (4.8). 
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For k* ■ (S«N)/2» thtt unlqua solution It glvon by y(x) ■ uhoro 
v(m) is givon by (4.S) and u(x) is tho liait of tha oonvargant saguanoa of 
functions 


( ^.90) u*(x) « h(x) ■ 


s-v«»> yp sinh k(l-x) ♦ a-v<<> y, ^j.nh Kx 
sink k 


( 9.9b) u*"*'(x) ■ h(x) ♦ 


t gH<x,o 

^0 


lk*u"(0 - F(|»u*(|))l d(> a • 0»l,***» 


uhara gnCx.^) it givan by (9.7b-l). 
Proof. Lot 


Ok; bliVvliHl. I'J 
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(9.10) 


E « max {a-v«o», a-^Mij 


Firstf ue shou that each member of the sequence ( 9.9) satisfies |u'”(x)| i EM* 
a 3 0 » 1 » •••, for all x « ( 0 » 1 1 . Observe that 


|u*(x)| i E 


lyol sinh k(l-x) ♦ |yi| sinh kx 
sinh k 


(9.11) 


sinh k(l**x) * sinh kx 

< E max {|yol»lyiD — 

sinh k 


cosh k((l/2)-x) 

< E n 

cosh(k/2) 


< E n, X € 10.11. 

Nou assume that u*(x) satisfies |uMx)| i EM for some i 2 0. 
Then from ( 9.9b). 


(9.12a) Iu«*nx)| < |h(x)| + ( gk(x.?) |k*u‘(f) - F(|.u»(n)l d?. 

^0 

Nou, if |u‘(x)l < EM. then |e''« « >u» (x) I « e''‘“» |u*(x)| < e''««>EM. Thus, if 
eV(Niyi(]() > 0, then 0 < f i (x.e''* “’u* (x)) i (Si*H|)e''‘**u‘(x) by hypothesis. 

This in turn implies from (9.5) that if q(x) i. 0. x c (0.11. then 0 ^ 

F(x.u‘(x)) < e-''‘>*»(Si*Ni)e''«>‘»u‘(x) + (& 2 >Nz)u‘(x) * (S*N)u>(x) ■ 2k<ui(x) if 
u‘(x) > 0. Similarly, it follous that 2k*u‘(x) = (&4’N)uUx) i F(x.u>(x)) < 0. if 
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u‘(m) i 0. Th»rttfor« th« tmt% |k*u>(|) - P((,uM|))| in thn Intngrnnd of (4.12n) 
will bo bounded by k*|u‘(()|» and it follows iron (4.12a) that 


|u>*nK)| i EM < 


tinh k(l>K) * sinh kn 
sinh k 


) * gk(K,|) dt k>EH 


(4.12b) 


sinh k(l-K) * sinh kx 

• Dt { ) ♦ (1 

sinh k 


sinh k(l'X) * sinh kx 
sinh k 


) Dt 


> Dt . 
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Thus* by induction* |u"*(x)| i CM* n > 0* 1* 2* •••. Now dafina 

(4.13) a"»*nx) » u»*nx) - u"»(x) 


and 


(4.14) lla"**'|| « max |a"**'(x)|. m ■ 0* 1* 

0 <x<l 

Since dfi/dy is continuous* dF/du is continuous so that wa can apply tha naan 
value theorem to F(x>u”*(x)) - F(x*u**'^(x)) to obtain 

f ’ dF 

(4.15) e^^nx) = \ 9k(x»5) Ik* (?*u"*(5)-d(5)e"'(5))) a"*(0 df* 

^9 du 


where 0 < d(^) < 1. Note that from (4.5) 


(4.14) dF/du * e*''«>‘» dfi/dy dy/du ♦ q(x) 


■ dfi/dy ♦ q(x). 

Also, if |e''*x’u(x)| < e''<*'*En* then 0 < Si S df i(x*a''< ***u(x))/dy i Ni* by hypothesis. 
This in turn implies from (4.16) that 0 £ Si-»St i dF<x*u(x))/du S Ni^Nt* if q(x) 2 0* 
or 0 ^ S < dF/du i N. Since ue have already shown that each mambar of tha sequence 
( 4.9) satisfies |u"*(x)| i EM* x e 10*1)* it follows that 

dF 

(4.17) 0 < S < — (5,u'"(f)-d(|)e"*({)) 1 H . 

du 
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Thui with th« oholc* 

(4. IS) 

th« braokttAd tar« in tha intagrand of (4.15) tatiaflaa |k* - sr/Su| i (N'S)/2 
and 


iiH 
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k> ■ (& « K) / 2 


1 t' 

|a<"«UK)| i - (N - S) \ 9k(x,$) d|- lla^ll . 

2 'A 


(4.19) 


1 ooah k((l/2)**K) 

• (l/2)(n - S) - (1 ) ||a"|| 

k* eoah(k/2) 


5 Hh . 


m > 0> 1» 


uhare 


(4.20) 


H - S 1 

flK » (1 ). 

N 4 S coih(k/2) 


Since (4.19) holds for all m c I0rl)» ||e'”^'|| S mh lla'"ll • Mote that 

Pk < 1> and ||e'"*^|| < (pk)"* lla^ll • Thus (u"*(x)} given by ( 4.9) is a Cauchy 

sequence in the space of continuous functions ulth the norm defined by (4.14). 

Uo can conclude that the sequence ( 4.9) converges uniformly co a limit u(x) 
which satisfies the integral equation (4.7a) and is such that |u(x)| i EH> 

K e to. 11. Then y(x) > e''*x'u(x) is a solution to (4.2a) with (4.2b-l) satisfying 
|ytx)l < e''‘*‘*Eh, X « 10, IJ. 

To establ.LSh uniqueness, let yi(x)> yz(x) be two solutions in which 
|yi(x)| < e''‘>‘’EM and |y2(x)| < e''’>“E«. x e 10. ll. Then ui(x) « e*''‘ •‘•yi(x) 
and uz(x) * e*''* ***yz(x) both satisfy the integral equation (4.7a) with 9k(x.() 
given by (4.7b-l)» h(x) given by (4.70'!) and k^ ■ (&4^N)/2 as specified in 
Theorem lA. By the same analysis that leads to (4.19), we can show that 
|ui (x)-uz(x) I < Pk llui-uzl| , X e 10,1). or 

llu, - Uzll i Pk llui - uzll . 


Since Pk is given by (4.20) and pk < 1. we must have ||ui'Uz|| • 0, 
or Ui(x) = Ui(x). Thus yi * y*. 

The importance of Theorem lA is that dfi/dy and f« are required to satisfy 
conditions only for all x c (0.1 1 and all y such that |y(x)| is bounded by an 
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•xprMsion dependant on boundary valuee and p(x). Of eouraa* the aKiatanoa and 

uniquanais applies only to functions that satisfy this oonstraint. This kind of 
s 

approach is uall notivatad though for problaM that arise naturally in the applied 
soianoas» sinoa y(n) usually raprasants a physioal quantity uhioh sMy be knoun to 
be bounded in absolute value or uhioh nay be of one sign. With this in nind the 
following additional result is obtained. It inposas conditions under the 
oonstraint that y(x) be nonnagativa and bounded above. In the interest of spaofi 
ua state it uithout proof. 

Thaoran IB. In the boundary value problem (4.2a) uith (4.2b>l) let ye» yi be 
nonnagativa and let max <yo» yi> S H. Suppose there exists an Ni > I and a 
6i 2 e such that 0 $ Si S Sfi/Sy $ Ni for all x c and all y such that 

0 S y(x) i x c 10,1 )» uhara E is defined in (4.10). Suppose further 

that 0 < fi(x,y) i (l/2)(Ni*Si )y for all x s 10,11 and all y such that 0 S y(x) 

£ a^"**En, X c (0,11. Than there exists a unique solution satisfying 
0 £ y(x) £ a''***’En, x c (0,1). Let S ■ Si^Sx and N ■ Ni«2Nx uhara N( and St are 
defined in (4.8). For k* ■ (S«N)/2, tho unique solution is given by y(x) ■ 
aV(xiu(x) uhara v(x) is given by (4.S) and u(x) is the limit of the oonvargant 
sequence of functions given by ( 4.9). 

The proof is similar to the proof of Thaoran lA and requires showing that 0 £ 
u*"(x) £ En. X s (0,11. m > 0, 1, •••. 

Similar theorems can be stated for (4.2a) together uith (4.2b-2) and (4.2a) 
together uith (4.2b-3i. Ua make note of one, for example, (4.2a) uith (4.2b-3). 
This problem can also be converted into the integral aquation (4.7a) uith 

1 f (k cosh k(l-S) * v^(l) sinh k(l-()) sinh kx, 0 £ x < |, 
(4.7b-3) gi<(K,() * - < 

K 1 sinh k( (k cosh k(l-x) * v''(l) sinh k(l-x)), ( < x £ 1, 


where 


K ■ k (k cosh k * v^(l) sinh k). 


and 

(4.7C-3) 


h(x) » e*''‘*» ye 


k cosh k(l-x) ♦ v''(l) sinh k(l-x) 
k cosh k * v^CD.sinh k 


Theorem 3B. In the problem (4.2a) together uith (4.2b-3) assume 0 £ ye £ H. 
Suppose there exists an Ni > 0 and a Si 2 0 such that 0 £ Si £ Sfi/8y £ Ni for 
all X € (0.1) and all y such that 0 £ y(x) £ x i -vt eiff, f (o,i|. Suppose 
further that 0 £ fi(x,y) £ (l/2)(Ni«Si )y for all x c 10,1) and all y such that 
0 £ y(x) £ e''* *’n, xe(0,l). Ifv^(l)2C, than there exists a unique 
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■olution tatisfying • $ y(N) S •viN)*vi«if|, n ^ |i»i|. dafifM S • 

•ltd N ■ Ni«2N| yh«T* Na and Sa ata fivan ky (4.d>. Than far k> • (N«ft)/2» tha 
unigua aolution is fivan by y(x> ■ a''*"*ti(N> iihara v(x) is givan by Cd.S) and 
u(N) is tha liait of tha eonvargant saguanoa of funotiona 

k eosh k(l**N) ♦ v^(l) aiidt kd-x) 

(4.21a) u*(x) • h(x) ■ a*''*** y^{ - - — — ' ) 

k OMh k * v^(l) sinh k 

(4.21b) u<"*'(x) ■ h(x) ♦ { gH(x,|) (k*ty»(|) - r(|.u*(|))l 4|. a ■ I, !#•••. 

uhara 0k(x,n is givan by (4.7b-5). — i« 

ORIGINAL la 

OF POOR QUALITY 

S. rUTURE ZNVCSTXOATXON 

It is of intak'«st to ac^an assuaptions for seaa of tha sosults alraady 
astablishad. For instanoa# ena oould oensidar raaoving tha assuaptions that 
q(x) of (4. Sc) ba nennagativa and tha assuaptions on fi. Thay oould ba ra- 
plaoad uith tha singla assuaption that thara axists a & 2 4 and an H > 

0 such that 

(5.1) 0 i 8 < ar/du i N, 0 i r(x,u) i (i/2)(n^8)u. 

for all u as constrained in Thaoraa 1B> for instance. Ua would also like to 
allow dfi/dy and/or dF/du to )>• nagativa, a.g.> to consider eases in 
which 

(5.2) -d S dfi/ay i N, d > •' 

Ua can state a simpla result for (4.2a) uith (4.2b-l) In iidtieh ua assuna that 
-d i dF/du < N for all ii, 0 i d ^ M/2» and taka k* • N/2. It will work if 
N < 8(cosh'*2)>. (tore general boundary conditions can also ba considered. 
However. Green’s functions and their analysis will ba aora conplioated. 

A generalisation to the problaa 

^*>.5) y^' ■ f(x,y,y^), 

in which we assume for instance that 0 S 8i S dfi/4y S Ni» and either 0 < St i 
df/dy' £ N| or -’N{ £ df/dy^ £ >82 £ 0> will bf^ harder. Ha propose to consider 
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(5.4) 


y'' - ki*y ± 2kty^ ■ f - ki*y t 2kiy^. 


ORIGINAL PAGE IG 
OF POOR QUALITY 


In this instanos a diffarant noni will hava to ba usad to astabllsh oonvarganoa 
of tha saquanca of suooatsiva ayproxiaations and in turn tha axistanoa and 
uniquanass, a.p.. aax (ki<|y(x)| * 2ktly^(x>|. Praliainary invMtlgatlon ravaals 
that unlika rasults rapo;*tad in Saetlon 3> oonvarganoa will dapand on tha 
aagnituda of Ni and Nj and oonvarganoa will not always ba guarantaad aaraly 
if af/dy and df/ay^ ara nonnagativa and boundad abova. 

Baeausa of axaapla 3 of Saotion (2.4)» and axaapla (A-4) of Bppandix k, it 
appaars that rasults for 


(5.5) 


/ yi'' \ /ii(yi) 

\ yz" J \<2(x»yi.yi.yi‘') 


Aiy(0) + Asy^(O) = Ci, Biy(l) ♦ BtV^d) ■ Cj, 


will ba useful. Unfortunately* a worm ooaplioatad integral aquation will have to 
be analired* namely one of tha form 


(5.6a) y(x) « H(x) * r6(x*() F(|*y(|)*y^(|)) d( 

T T 

where H > (h^* h{) , F > (fi* fj) > and 


(5.6b) 


/ gii(x,|) gi2(x*t) 
6(x*H * ' 

\ g2i(x»5) g22(x.f) 

\ 


6. sons APPLICATIONS AND S0T1E NUI1ERICAL RESULTS 

Lets demonstrate an application of tha rasults reported in Section 3 to* for 
example* problem (?.l), i.e.* 

(6.1a) y" * KRy" - xRy', 

(6.1b) y(0) s 0* y(l) « 1. 

For this specific case of (4.2a) with (4.2b-l)* we hava 


la 


( 6 . 2 ) 


fi(x*y) B KRy" 


p(x) ■ -xR. 


OWGtNAL PAGS 13 

In tsriMi of (4.4) and (4.5), (4.1) baowMs OF QUALITY 

(4.3a) u"(x) ■ F(m,u(x)) ■ ••''‘•‘•KR(a''*>‘*u(M))'* ♦ q(x)u(x), 

(4.3b) u(0) ■ 0, u(l) ■ ■ axp(S/4), 

uhara 

(4.3c) v(x) ■ -x*R/4, q(x) ■ (-xR/2)* ♦ R/R, U(x) ■ a’''***y(x). 

For k* # 0, (4.3) can ba eonvartad into 

sinh kx 

(4.4a) u(x) n + \ gk(x.p<k<u(() - F(|,u(?))> d?. 

sinh k 


uhere 


1 I sinh k(l-() sinh kx, 0 < x < (, 

(4.4b) gk(x.5) S 

k sinh k I sinh k(l-x) sinh k(, ( < x i 1, 

In this problem, one expects that 0 S y(x) i 1. Ha shall daaionstrata 
a direct application of Theorem IB, and shou constructively that there exists a 
unique solution satisfying 0 < y(x) i e''* ^ » » exp(R(l>x*)/4), x 6 (0,11. 

The conditions of Theorem IB are satisfied as follous. 

(a) For this problem, fi(x,y) = KRyn, so 

(6.5a) dfi/dy = nKRy"*^ 

and 

(6.5b) 0 < dfi/dy < nKR(exp(R/4))"*l 

for all x e (0,11 and all y such that 0 < y(x) i ^ c (0,11. Note that 

(6.4a) * 0, Ni < nKR(exp(R/4))"*i 

in this case, while 

(6.6b) Sz = min q(x) * R/2, Nz * max q(x) * R</4 * R/2. 

0<x<l 0<x<l 

As defined in Theorem IB, 
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(ft.(o) 6 - R/2. N - nKfUaMpCRM))***' * 2(R>/4 *■ S/2). 

(b) If M < t0»ll» and y is suoh that 0 $ y(x) i n ^ than 

(«.7) 0 i fi(K»y) i (n/2) KR(aMp(Il/4))n-iy • (i/2)(Ni*Si )y» n 2 2. 

hceerding to tha eonelusion of Thaoraa 1B> for 
(6. ft) k> - (l/2)(8 ♦ N) > (l/2HnKR(axp(S/ft)n*l ♦ RV2 * (S/2)P) 

the saquenca 
(6.9a) 

axp(R/4) sinh kx 

C6.9b) u"*Mx) » ♦ \ gk(x.?){k*u"*(0 - r(C,u«*(|))}d(^. 

sinh kx 

as 0» has a limit u(x)» and y(x) ■ a''*’<)u(x) • axp(*Rx*/6)u(x) is tha 

uniqua solution of (6.1) satisfying 0 < y(x) i axp(R(l-x‘)/4)» x < I0»1I. In 
particular for tha saquanea (6.9)» 


Q) 


f- 


u*(x) > axp(R/6) sinh kx / sinh k» 


OF 


»■« 




||uin«1 - um|| < Hum _ u***^|| > 

uhere 

nKR(exp(R/6 ))•*-’ ♦ R*/2 ♦ R/2 / 1 \ 

(6.10) ilk ■ 1 

nKR(exp(R/6))»'-< ♦ R*/2 ♦ (3/2)R ^ cosh(k/2) / 

To obtain a numerical solution of tha limit u(x) of tha saquanea (6.9)» it 
can be approximated by the discrete solution Uo « u(0), ui * u(xi>> 
up.i s u(xp.i)> up * u(l) on a uniform grid* h ■ 1/p* Xi • h* xj ■ 2h* ***» 

Mp.i s (p-l)h* Kp s 1. By using the trapezoidal rule on tha integrand* U| can 
be computed by the expression 


(6.11) U| s h(X| ) 


•? 


0fi gk(Xi*xi) lk*Ui 


j=o 

uhere Op ^ op > h/2, Oj = h* j = 1* •••*?-!. 
should yield accuracy on the order of h*. 


- F(xi*U|))I* i * 0*1*>‘** p* 
Keller ( 1 ) has shoun that this 
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To eonputo tho approximations a saquanea of nat functions 

m > 9 . 1 > , can ba daf inad as 

OF POOR 

(6.12a) Ui* ■ h(Xi)f 

(6.12b) ■ h(X|) * ^ 9k(*i»Kj) Ik*uj* - F(X|»Uj*)|> i ■ 9>lf***> p> 

j»0 

uhare Oo ~ orp ■> h/2. tt| « h. j « It ***t p-1. By argunants similar to thosa 
given by Keller (II for (3.1)# one can shou that the limit of (6.12) exists and 
is the unique solution of (6.11). Also# as p the contractive paraMter 

for the sequence (6.12) will converge to the contractive parameter t>k for the 
corresponding continuous sequence. 

The following table indicates some values of the contraction parameter Pk 


QUAuny 


n 

R 

K 

Pk (6.10) 

2 

10 

.01 

.76 

2 

10 

.1 

.S2 

5 

10 

.01 

,999 

5 

10 

.1 

.9999 


Table 1. 


It is important to realize that the value of Pk is merely a bound on the 
contraction of the sequence. In practice, the actual number of iterations 
required to achieve a given error tolerance way be significantly less than the 
number expected by the value of pu- To illustrate this ue programmed (6.12) for 
problem (6.1). Programming was done in standard Fortran on and IBI1 3033. For 
example consider the case n=2# R=10# and K=0.01 in Table 1. Assume that p is 
sufficiently large that the contractive parameter for (6.12) is nearly Pk of (6.10). 
Then if ||u’ - u®|| < 1# the relation Hu"*^ - u"*|| < (Pk)"* llu’ - u®|| 
implies it should take at most 33 iterations to get ||u"**^ - u"*|| i 10*'. It 
actually took m^l = 16 with h - (1/50). Consider the case n^S# R«10# and K^Q.D!. 
Note that pn is very close to 1 in this case. In such a case# the rate of 
convergence of the sequence may be very slow# thus causing the number of iterations 
required to achieve a given error tolerance in the programming of (6.12) to be 
large. Keller (II suggests Newton's method as one alternative to a more rapidly 
convergent scheme. However, analysis such as the following can also be applied. 


16 


First* ue obsarvs that tha initial approKimation u*(x) of (*.9a) is ona 
suggastad by tha intagral aquation itsalf. If tha initial approxination is 
chosen closer to the limit of tha saquanca* tha number of iterati<xts required to 
achieve a given error tolerance will be reduced since ||u^ - u^|| mill be smaller. 

To this end ue can appeal to comparison theorems C2I. For instance* if yt(x) is 
the solution for n«2» R«10> Ks.Ol* and 0 i ya^K) i 1* then ya(x) i ys(x) uhera ya is 
the solution for n«5* R«10* K>.01. Thus in the case n>S* R«10* K>.01* the initial 
app.roximation in the sequence (6.9) may be taken as u*(x) ■ exp*''<>'>y 2 (x} instead 
of u*(x) given by (6.9a). 

Secondly* ue suspect that the solution y(x) « e''<x>u(x) * exp(-x*R/6) u(x) 
satisfies 0 < y(x) i 1 for problem (6.1). If ue can shmi that 0 < exp(-x^R/6)uai(x) 

< 1 for each m* then 

dF 

(6.13) R/2 < — (K*u"*) < nRK ♦ R*/6 ♦ R/2 

du 

for each m. Thus if ue choose 

(6.16) k* « (1/2) (nRK ♦ R*/2 ♦ R) 

the sequence (3.33) uill contract according to 

||U»*1 - u»|| < fiH llu" - u"*-'|| 


OF POOR QUALITY 


uhere 


nRK + R*/2 1 

(6.15) = (1 ) 

nRK + R*/2 + H cosh(k/2) 


a..J k is given by (6.16). Compare the following Table 2 to Table 1. 


n 

R 

K 

iik (6.15) 

5 

10 

.01 

•73 

5 

10 

.1 

.75 


Table 2. 

When we us<>d k^ given by (6.16) and the solution for the case n«2* RslO* K« 
.01 as the initial approximation in '*'he sequence (6.12)* it only took 6 iterations 
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uith h ■ 1/SO to got i|u*** - u*|| i 10*^ for tho ooso Salt* in ptobloM 
(0.1). For tho COSO naO» Ral0» Ka.Olt uo usod k* glvon by (0.14) and tho 
s.»lution for tho oiko n«S» R>10» K*.01 os tho Initial opproxiootion and it only 
took 3 itorations. Finallyf for tha oaso naS> Ral0> K*.01> no usod k* given 
by (0.14) and tho solution for tho oaso n«4, Ral0> Ka.Ol os tho initial 
approxiaation* and it only took 3 itorations. 

Considor tho voetor probloa (2.4) shown in Sootion 2> 


original F.^ V 

OF POOR Q- >•' 


(O.lOa) 



■ f”’ 

S(xyj-' ♦ lyt - yiyi') 



(O.lOb) 


yi(0) ' 0» yj(0) a 0» 


yid) a 1> yi'(l) a 0. 


There are several ways to extend tho approach to systoas. For instanco* uo noto 
in this case that 


(4.17a) dfi/dyi ■ 0, dfi/dy^ « 1, dfi/dyi' » 0, dfi/dyj'’ ■ 0. 

(4.17b) dfx/dyi a -sya'* afj/dyi a 3s, dfi/dy,-' a o, dfa/ay*' a s(x-yi)» 


Then* write (4.14) in terms of ui = yi and ut * e'*»’‘y 2 * 

(4.1«a) D*ui - li*U 2 a -Fi 

(4.18b) -mi*ui + (D* - n*)u 2 » -F 2 . n* » ni<-»n 2 ‘ 

where D* * d*/dx* and 

(4.18c) Fl a ll*U 2 - e-*'*>‘U 2 . 

(4.18d) Ft a mi*ui ♦ (ni*U2 “ 2n2(U2'-n2U2) - e'****f2). 

The boundary conditions are 

(4.18e) Ui(0) a 0> U2(0) a 0* Ui(l) a 1, Ui^(l) a Q. 

Problem (4.18) can be converted into an integral equation of tho form 


(6.19) 


I Ui(k) \ 
\ U2(x) / 

V ' 


' hi (x) 
\ h2(x) 



G(x*u 


Fi(?.U2(5)) 

r2(5.uid)*u,(n»u,'(?)) 


dr 


For explicit representation of hi* h 2 * Fi* F 2 * and 6 soo (8). 

To get a numerical solution to (4.19)* one can try programming tho soquonoo 
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(6.20«) ui*(x) ■ hi(x)f ■ hi^(x)* 

(6.20b) ut*(x) • hz(x)> Ui*'(x) • hs^(x)t 


ORiaiNM- P.Vf a 

Of fOOR QUALITY 


(6.20c) 


'Ui‘-*'(x) 

luj*^'''(x) 



\ /rid.UiMD) \ 

♦ \ G(x,5) df, 

/ ^0 \r2(^Ui‘'(|),U,^|),U2»•'(5))/ 



i: 


Fl(|rU2^(()) 

6^(X>|)( /d|. 

‘r2(|.UiM0.U2^(0»U2^'(|))/ 

1 ■ 0 » 1 » • • • . 


Suppose ue take S > 0» and try to find a unique solution for uhioh 
(6.21) 0 < yi S 1, 0 < yi' S P, -Q S Vz S 0. -M S Yt' i 0, 

for all X € 1 0 > 1 1 . Ue propose choosing values 

(6.22a) ni* » 5S/2> 2nt • max| (x-yi )S|/ 2, 

(6.22b) mi* * max (-Syz'') / 2. li* » (I ♦ a*"») / 2 

Then* if each member of the sequence (6.20) satisfies the bounds (6.21)» and if 
the sequences converge uniformly uith respect to some norm* ue can conclude 
that a unique solution to (6.16) satisfying the bounds (6.21) exists. Its the limit 
of the sequences (6.20). This has not been done rigorously yet. Houever* ue 
shall report some numerical results uhich are encouraging. 

The sequence (6.20) uas replaced by an approximation in uhioh ue evaluated 
ui®, U 2 ®. and U 2 ®' discretely at p points X| = (j-l)h. j « l.***.p» 
h - l/(p-l). Then ue obtained the values of the remaining members of the 
sequences at each discrete point by evaluating the integral using the 
trapezoidal rule. The interation uas stoped if and uhen 


max { |ei"»(X| )| » I (ei"*(xj ) )■' I * |ez"*(xj ) I * I (ez*"(xi ))'' |) < lO'^ 
3 

where 


ei"*(X| ) » yi"*(xj )-yj"*'Uxj )* (ei'"(X|))'' * (yi"»(xj )-yi"**'(xi))'» i ■ 1*2. 


The following results were computed using standard Fortran on an IBM 370 3033. 
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eas« 1 ou« 2 


oMa 3 


s 

.11576 

2.081 

6.868 


ni 

.4167 

1.767 

3.210 

OWQINAL 


.006682 

.2601 

.5151 

OF POOR QUALITY^ 

«i 

.4167 

2.04 

6.934 



.9978 

.9410 

.8937 


h 

1/50 

1/50 

( 

1/40 


iterations 

4 

6 

11 


CPU time 

46sec 

77sec 

93sac 


yz(l) (computed) 

-3.0622 

-3.968 

-5.700 


yzCl) value 
reported in (5) 

-3.0622 

-3.961 

-5.503 



C0NCLUDIH6 RENARKS 

Analytical Advantagaa and Disadvantages . 

Constructive existence and uniqueness is given for a large class of iion> 
linear problems. Houever* explicit Green's functions can be difficult to con- 
struct and analize. Also« analysis for systems can be complicated. 

Numerical Advantages and Disadvantag es. 

One advantage is that boundary values are built right into the integral 
equation. In numerical procedures > one does not necessarily need to approximate 
the derivative at one end if only the unknown is given as in initial value 
methods. No derivatives need to be approximated and therefore no difference 
approximations are required. Ue know before a numerical nwthod of approximating 
the sequence of successive approximations is applied that the sequence converges 
to the unique solution of the problem. 

One disadvantage is that in practice a large number of mesh points nay be 
required to approximate the function sequences accurately. Another is that slow 
convergence may be a problem. 
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AppandiM A 


In a problam oonoarning radiation haat tranafar for annular fins lfl» tha 
tamparatura distribution is shown to ba govarnad by tha anargy aguation writtan as 


(A-la) 


y"’(x) ■ q(x) (y(x))' ♦ p(x) y'(x), 0 S x £ 1» 


(A-lb) 

where 

(A-lc) 

and 

(A-ld) 


y(0) « 1.0. 


y'(l) » 0. 


g(x) 


0 

(l-x) tan a * $ 


p(x) 


tan tf 1 

(l-x) tan Of ♦ d X ♦ p 


The quantities or. d. and p are constants related to angle of taper, fin 
thickness at tip. and radius of fin base and fin tip. 

A flat plate nodel of a catalytic converter (10) leads to tha nonlinear two 
point problem 

(y')» 

(A-2a) y" , 

(2 - y) 

(A-2b) y(0) = yo. y(l) = 0. 

The unknown y represents mole fraction of a gas. 

A problem involving fluid flow in a two-dimensional channel (11) leads to 
the following boundary value problem governing the velocity distribution 


(A-3a) 


f"" = R(ff"' - f'f^') 


(A-3b) f(0) » 0. f"(0) » 0. f(l) « 1. f'd) « 0. 

Me can write this problem in vector form with yi > f and yj ■ f", 

yi 

R(yiy?' - yi'yz) 



01 


(A-«b) 


yi(0) * 0. yz(0) * 0 


yi(l) ® 1. yi^(l) ■ 0. 
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